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Abstract: Based upon flavor SU(3) symmetry, we perform global fits to charmless 
B decays into one pseudoscalar meson and one vector meson in the final states. We 
consider different symmetry breaking schemes and find that the one implied by naive 
factorization is slightly favored over the exact symmetry case. The (p, ff) vertex of 
the unitarity triangle (UT) constrained by our fits is consistent with other methods 
within errors. We have found large color-suppressed, electroweak penguin and singlet 
penguin amplitudes when the spectator quark ends up in the final-state vector meson. 
Nontrivial relative strong phases are also required to explain the data. The best-fit 
parameters are used to compute branching ratio and CP asymmetry observables in 
all of the decay modes, particularly those in the decays to be measured at the 
Tevatron and LHC experiments. 
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1. Introduction 

Thanks to the B-factories, a plethora of data on rare hadronic B meson decays have 
become available in recent years. Because they involve VT-mediated charged-currents 
through mixing and/or decay, these decay modes provide particular useful informa- 
tion on the CP-violating weak phases and magnitudes of elements in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [|l], |^ for the quark sector of the standard model 
(SM). Advances in both experiment and theory have helped us narrow down these 
parameters to a high precision. Through such efforts, it therefore becomes possible 
for us to search for evidences of new physics, if any. 

Due to the hadronic nature of particles involved in the decays, strong phases 
associated with the decay amplitudes that are derived from short-distance physics 
as well as final-state interactions are also important. Even though they cannot 
be computed from first principles, these phases play a crucial role in direct CP 
asymmetries. Determination of their pattern and magnitudes in B decays give a test 
to our knowledge of strong dynamics in the SM. 

An approach utilizing the flavor symmetry to relate magnitudes and strong 
phases of amplitudes 0, ^, ^, H, has been taken to analyze the rare B decay data. 
It has the advantage of reducing model dependence for computing matrix elements 
of hadronic transitions, in comparison with the usual perturbative approaches. 

In Ref. IQ, we have updated the analysis for B decays into two charmless pseu- 
doscalar mesons in the final states, and further tested the flavor symmetry assump- 
tion by considering several different breaking schemes in the amplitudes. By perform- 
ing global fits, we find that our results are robust against fluctuations of individual 
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data with large uncertainties, and different schemes have roughly the same predic- 
tions. 

In this article, we concentrate on the rare B VP decays, where V and P 
denote charmless vector and pseudoscalar mesons, respectively. There have been 
some numerical works in the perturbation framework of quantum chromodynamics 
(QCD) to calculate the decay rates and CP asymmetries of these decays over the 
years. Naive and general factorization analyses were considered in Refs. |P, 0, |ll 



The QCD factorization (QCDF) method was employed in Refs. |T|, |T3|, y, [T|, |T|, y]. 



The calculations using the perturbative QCD approach are scattered in Refs. [|18| [19 



20| , |2T| , p2| , p3l . Recently, the soft-coUinear effective theory (SCET) was also used 



in Ref. [^. In parallel, some attempts that apply the flavor symmetry to the VP 



decays are given in Refs. [^, |2^, |2^, ^ . 

The B VP decay modes present a richer structure than the PP final states 
because the light spectator quark in B meson can end up in a spin-0 or spin-1 meson, 
even though the quark-level subprocess is exactly the same. Moreover, the number 
and precision of observables in these modes (particularly the strangeness-changing 
ones) have improved considerably in recent years. Totally, there are 52 observables 
in the VP decays. All the branching ratio and CP asymmetry observables in the 
strangeness-changing decays of B^'~^ mesons have been measured. The branching 
ratio of p~^K^, in particular, provides valuable information on the magnitude of one 
type of QCD penguin amplitude. In contrast, the observables in the strangeness- 
conserving transitions are mostly measured in the B~^ decays. Moreover, some data 
points have shifted by noticeable amounts. For example, the central values of the 
branching ratios of B~^ — > p~^i]^'\ B^ K*~^r], and K*'^7i~ have dropped by about 
30% from five years ago. The branching ratios of B^ — > p^n^ also move significantly 
upward and downward, respectively. Therefore, we consider it timely to re-analyze 
the data and, at the same time, relax some of the assumptions made in Ref. p8|] in 
view of the better data pool, and make predictions for the Bg decay modes which 
are going to be measured at Tevatron and LHCb. 

The structure of this paper is as follows. In section ^ we introduce the notation 
used in our approach and present both measured observables and amplitude decom- 
position for the decay modes. In section ^, we show our fitting results of the theory 
parameters in different schemes. Discussions and predictions based on our fits are 
given in section |[ Section ^ summarizes our findings in this work. 

2. Formalism and Notation 

For a two-body B —>■ V P decay process, the magnitude of its invariant decay am- 
plitude M is related to the partial width in the following way: 
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where p is the 3-momentum of the final state particles in the rest frame of the 
B meson of mass m^. To relate partial widths to branching ratios, we use the 
world-average lifetimes r+ = (1.638 ± 0.011) ps, r° = (1.530 ± 0.009) ps, and Tg = 
(1.437 ± 0.031) ps computed by the Heavy Flavor Averaging Group (HFAG) [p9|] . 
Each branching ratio quoted in this paper has been CP- averaged. 

To perform the fiavor amplitude decomposition, we use the following quark con- 
tent and phase conventions for mesons: 

• Bottom mesons: B^ = db, = bd, B^ = ub, B^ = —bu, Bg = sb, Bg = bs; 

• Pseudoscalar mesons: 7r+ = ud, 7r° = (dd — uu) / -\/2, vr^ = —du, = us, 

= ds, K = sd, K~ = —su, i] = {ss — uu — dd)/^/3, v( = {uu + dd+2ss) /^/6; 

• Vector mesons: p"*" = ud, p° = (dd — uu)/^/2, p~ = —du, lu = (uu + dd)/V2, 
K*+ = us, = ds, K*^ = sd, K*- = -su, = ss. 

The 7] and rj' mesons correspond to octet-singlet mixtures 

= ?78 cos 6*0 - r^i sin 6^0 , (2.2) 
77' = sin 6^0 + ?7i cos 6^0 . (2.3) 



As shown in Ref . ||2^ , varying the mixing angle Oq does not improve the quality of fits. 
For convenience, we fix 9q = sin^^(l/3) ~ 19.5° according to the above-mentioned 
quark contents of 77 and rj'. 

We list fiavor amplitude decompositions and averaged experimental data for 
B —>■ VP decays in Tables |I| and 0. Values of measured observables are obtained 
from the latest 2008 summer results of the HFAG |p9| . 

In the present approximation, we consider only five dominant types of inde- 
pendent amplitudes: a "tree" contribution T; a "color-suppressed" contribution C; 
a "QCD penguin" contribution P; a "fiavor- singlet" contribution S, and an "elec- 
troweak (EW) penguin" contribution Pew- The first four types are considered as 
the leading-order amplitudes, while the last one is higher-order in weak interactions. 
Depending upon which final state meson the spectator quark in the B meson ends 
up in, we further associate a subscript P or y to the above-mentioned amplitudes. 
For example, Tp and Ty denote a tree amplitude with the spectator quark of the B 
meson going into the pseudoscalar and vector meson in the final state, respectively. 
These two kinds of amplitudes are different in general. In the following, we will 
suppress the subscripts P, V when discussions apply to both classes of amplitudes of 
each type. 

There are also other types of amplitudes, such as the "color-suppressed EW pen- 
guin" diagram P^w^ "exchange" diagram E, "annihilation" diagram A, and "pen- 
guin annihilation" diagram PA. Due to dynamical suppression, these amplitudes are 
ignored in the analysis. 
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Mode 


Flavor Amplitude 


BR (xlO^^) 


Acp 




K 


PP 


0.68 lb 0.19 








Pv 








p TT 


~^i'^v + cp + Pv - pp) 


o 7+1.0 


n 07+0-12 




p TT 


— :^(*p + Cv +pp — Pv) 


in 

J^U-9_i 5 


n 09 -1- n 1 1 

U.UZ zt U. 11 




P+T] 


~V5^ + cy +PP +py + Sv) 


n 1 in 


nil 1 nil 
U. 11 =t U.ll 




+ / 

P V 


1 / -L 1 1 1 \ A \ 

-^{tp + CV +PP+PV + 4sy) 


n 1 +3.7 
9-1-2.8 


—0.04 ± 0.28 






-^{W + cp + PP + Pv + 2sp) 


6.9 ± 0.5 


-0.04 ± 0.06 




1 4- 


Sp 


< 0.24 








PP 




- 




K*^K 


Pv 


< 1.9 








-{tv+Pv) 


16.42 ± 1.96'=' 


0.12 ±0.06" 
-0.04 ±0.13'' 






-{tp+pp) 


7.58± 1.25'^ 


-0.14 ±0.12"^ 
0.06 ±0.13" 







1 / 1 N 

~2\^P +CV - pp -pv) 


2.0 ± 0.5 






P^T] 


■^{cp - Cv - PP - Pv - Sv) 


< 1.5 






/ 


1 / A \ 

-^(cp - CV - PP - Pv - 4sy) 


< 1.3 






n 


2(Cp - Cv + PP + Pv + 2sp) 


< 0.5 






UJTj 


— -^(Cp + Cv + PP + Pv + '^Sp + Sv) 


< i.D 






UJff 


^^(cp + Cv + Pp + Pv + '^sp + 4sv) 


< i.y 






(pTT 




^ n 98 

U.ZO 






4>v 




< 0.52 






4>v' 




< 0.5 




Bs - 


p-K+ 
p^' 

ir% 
ir\' 


-j^icv-pv) 

-{W +Pv) 
-{tp+pp) 

-^(cp-pp) 
— \=(cv — VP + Pv + Sv) 

^ ( r'-^r -\- Ot) d -\- TlT / -1- 4- ^ 

^\ y ' ^PP pv 1 ^'^V ) 
-^{cp+pp + 2sp) 




- 




0" 


Pv + sp 







" Values obtained using the method described in Ref. [p8| 



Table 1: Flavor amplitude decomposition and measured observables |]30| , 31, 32, ^ of 
strangeness-conserving B — > l/P decays. The time-dependent CP asymmetries A and S, 
if applicable, are listed in the first and second rows, respectively. 



The QCD penguin amplitude contains three components (apart from the CKM 
factors): Pt, Pc, and P„, with the subscript denoting which quark is running in the 
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Mode 


Flavor Amplitude 


BR (xlO"^) 


Acp 


rS^ - 




Pp 


1 n ri 1 no 


-0.020_oo61 
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4-9-1.9 
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6.7 ± 0.5 


0.02 ± 0.05 




(t)K+ 


Pp + Sp 


8.30 ± 0.65 


0.034 ± 0.044 


- 




-{t'p+p'p) 


10.3 ± 1.1 


-0.25 ± 0.11 






^,ic'y-p'p) 


2.4 ±0.7 


-0.15 ±0.12 




p~K+ 


-{i'v+Pv) 


8.6l";^ 
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p J\ 


-T^^'^p-Pv) 


t; /i+o!9 

5-4_i.o 


n no _L n on 
— U.Uz zb u.zy 
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j\ rj 
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1 c; q _|_ 1 n 


vj. xy ziz u.Utj 






■^{dy+Pp + 2py+ASy) 


3.8 lb 1.2 


—0.08 ± 0.25 






^(dp+p'y + 2s'p) 


5.0 ±0.6 


0.32 ±0.17 








0.45 ± 0.24 






p'p + s'p 


8.3^1:^ 


0.23 ±0.15 

^-"^"^-0.18 


Bs- 


K*+K- 


-ip'p + t'p) 
-(Pv + iv) 
p'p 

rl 

Pv 








p rj 


v& p 

1 r' 

v^^p 














UJTj 


y/^y p ^ ■^^ p> 








i < ^tJ 

UJ 1 1 


]-(c' + 2s' "1 




















-^{p'p + P'v - 4 + «p - 4) 










-j^{2p'p + + cV + 2s'p + As'y) 







Table 2: Flavor amplitude decomposition and measured observables pO| , 31, 32, ^ of 
strangeness-changing B VP decays. The time-dependent CP asymmetries A and S, if 
applicable, are listed in the first and second rows, respectively. 



loop. After imposing the unitarity condition, we can remove the explicit t-quark 
dependence and are left with two components: Ptc = Pt — Pc and Pm = Pt — Pu- For 
simplicity, we assume the t-penguin dominance, so that Ptc = Ptu = P- The same 
comment applies to the EW penguin and singlet penguin amplitudes, too. 

In physical processes, the above-mentioned flavor amplitudes always appear in 
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specific combinations. To simplify tlie notations, we tfierefore define tlie following 
unprimed and primed symbols for AS = and \AS\ = 1 transitions, respectively: 

J. VUrp jyu I yC \ pC J./ T (V^ _L \ 

^ — -'dft-^ K^db "T ^dbJ-'^EW ' — -'sfeSi-i yi "r J- sbJ-'^EW ? 

c = YdbC - ^db + , c' = Y^bicC - {Y^b + Ygb)PEw , 

p = -(17, + Y^b) (P - ^^fiy) , P = -(^6 + Ysl) {^^pP - IpEw^ , 

s ^ -{Y^b + Y^b) [s - \pew) , s' ^ -{Y^b + Y^b) [isS - \pew) , (2.4) 

where Y^^^ = VqiqV*,b {q G {d,s} and g' G {u,c}). Here we also keep the P^w 
amphtude for completeness, though it is ignored in the subsequent analysis. Again, 
all the above amplitudes are to be associated with subscript P or V, depending on 
the process. Here we have explicitly factored out the CKM factors, but leave strong 
phases inside the amplitudes. 

From AS = to \AS\ = 1 transitions, we put in SU(3) breaking factors 
iTpy.icpy^ and ipp y for Tpy, Cpy, and Ppy, respectively. If some type of ampli- 
tudes is factorizable, the corresponding SU(3) breaking factor is either fx/ f-K = 1.22 
or fK*/fp = 1.00 p3. For example, we have for the —>■ K*^7i^ decay: 

AiK*+n-) = -Y^b^tTp + (Y^b + Y:b) i,Pp . 

This can be obtained from the complete set of flavor amplitude decomposition given 
in Table 0, Table and appropriate forms of Eqs. (pl^). 

In this analysis, the CKM factors are expressed in terms of the Wolfenstein 
parameterization to O(A^). Since A has been determined from kaon decays to a 
high accuracy, we will use the central value 0.2272 quoted by the CKMfltter group 
^ as a theory input, and leave A, p = p(l — A^/2), and f] = ri{l — A^/2) as fltting 
parameters to be determined by data. 

For the B meson decaying into a CP eigenstate /cp, the time-dependent CP 
asymmetry is written as 

j,^^^,^-nB^-fcp)-V[B^^fcp) 



T{B^^ fcp) + T{B^^ fcp) 
= S sm{AmB ■ t) + A cos{AmB ■ t) , (2.5) 

where ArriB is the mass difference between the two mass eigenstates of B mesons 
and t is the decay time measured from the tagged B meson. 

3. Fitting Analysis 

In this section, we present the following two schemes in our flts: 
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1. exact flavor symmetry for all amplitudes {i.e., C,Tpv — ^Cpv ~ ^Ppv — l)j 

2. imposing partial SU(3)-breaking factors on T and C amplitudes only {i.e., 

iTp,Cp = fK*/fp and ^Tv,Cv = /i^//7r, while ^Pp,Py = 1); 

We have assumed exact flavor symmetry for the strong phases to reduce independent 
parameters in our fits. Besides, Tp is fixed to be real and positive in our phase 
convention {i.e., 6tp = 0). All the other strong phases are measured with respect to 
it. 

We further divide our fits into two classes: (A) the VP modes that do not involve 
singlet penguin contributions, and (B) all of the VP modes. As shown in Table |l] 
and Table 0, the modes that contain the singlet penguin amplitudes are those having 
rj, rj' , (j), or u in the final states. 

It is appropriate to list some major differences between the current analysis and 



Ref . p8[ . Throughout this analysis, we do not assume any strong phase relation be- 
tween the EW penguin, singlet penguin, and the QCD penguin amplitudes. Neither 
do we assume any strong phase relation between the color-suppressed amplitudes 
and the tree amplitudes. The relative size and phase of Pp and Py are always kept 
free. Moreover, we do not assume Sp to be small enough for omission. Instead, we 
keep and constrain its magnitude and phase. 

In the following, we perform fits to the observables in the B VP modes 
as well as \Vub\ = (4.26 ± 0.36) x IQ-^ and \Vcb\ = (41.63 ± 0.65) x IQ-^ for the 
above-mentioned two schemes. The inclusion of \Vub\ and \ Vcb\ helps fixing the values 
of A and \/p'^ + ff. However, we drop the branching ratio and direct CP asymmetry 
of the — i> K*^tt'^ decay from the fits because currently the BABAR Collaboration 
and the Belle Collaboration have a large disagreement in the branching ratio, whose 
weighted average is (2.42 ± 1.16) x 10^^ with a scale factor S = 1.77. As we will see 
later, our predictions based on best fits deviate much from these two observables. 

The fit results of theory parameters are summarized in Table |^. As given in the 
table. Scheme 1 of exact SU(3) symmetry is slightly worse than Scheme 2. As defined 
above, the main difference between these two schemes is in the scaling behavior of 
Tp and Cp between the strangeness-conserving and strangeness-changing modes. 
We have also tried other schemes, such as having additional symmetry breaking for 
amplitude sizes. However, either the fitting quality becomes worse or they involve too 
large SU(3) breaking (over 30%). We will present our plots and predictions mainly 
for Scheme 2. 

Some general features are observed in these fits. The two types of tree amplitude 
have roughly the same strong phases, with Ty larger than Tp by about 50%, largely 
driven by the branching ratios of p^vr^. The Cy amplitude is 3 to 7 times larger 
than the Cp amplitude. Both of them have sizeable strong phases relative to the tree 
amplitudes. Moreover, the strong phases of Cp changes abruptly when we enlarge 
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Parameter 




Sdieme 






1 A 


2A 


IB 


ZD 


Tp 


0.721 ± 0.088 


n TOT _L n non 

0. /z/ ± u.uey 


n Totx _L n nno 

0. / So ± o.oys 


0.791 ± 0.100 


Tv 






-1 -1 /^o4-0 131 


1 1 '7n+0.133 
l-170lo.ll8 


Stv 


1.9 ± 5.7 


2.3 ± 5.7 


0.4 ± 5.5 


0.6 ± 5.4 


Cp 


0.0931°:!^ 


n 1 o /I _u n OOQ 
U. io4 ± (j.ZZo 


n -r 70+0.138 
''-'-0.108 


^ ^ r-vrA-l-O 1 9"^ 
U.±ZZ_Q Qgg 


OCp 


lion _i_ A 


—107.7 ± 31.0 


133.0 ± 34.2 


1 /tn n-l-72.4 
149.0_43.o 


Cv 


^ ,^oo+n 226 

0.688l^:l7^ 


^ /(+0 209 

0.624l^:f^^ 


0.945 ± 0.142 


0.892 ± 0.139 


Scv 


-66.0^^^:^ 


-57.01^^:^ 


of-> n + 12 

-82.oii^:° 


-75.91}^:? 


Pp 


0.084 ± 0.003 


0.084 ± 0.003 


0.085 ± 0.003 


0.085 ± 0.003 


Spp 


-3.9 ± 10.2 


-5.7 ± 10.0 


-1.0 ±8.0 


-2.6 ±7.8 


Pv 


0.065 ± 0.004 


U.UDo zt U.UU4 


U.UDo ± U.UU4 


0.066 ± 0.004 


OPy 


i / i. / ± o.i 


172.6 ± 7.7 


172.2 ± 7.1 


1 TO K _i_ fi n 
i / Z.O ± 0.9 


Pew,p 

OPew,p 


„ „„„-i-n nnq 
0-0391°:°?? 


U.UOJ_Q Q]^Q 

1 +10.4 

■J '-'•-'--11. 9 


fin n+10.0 


„ „„ -1 +0 01 

0-0311°:°}° 

crn n+10.5 


7") 

^Ewy 


U.UoT ± U.U49 


0-052_oo4i 


n nnfi-t-0.027 
0-096_o.o30 


A r\0'~7 1 A AAA 

0.087 ± 0.029 


"Pew,v 


no ^7+52.0 


nn r>+82 

-90.21°^:° 


-113.5l|^ 


1 1 1 n+10.4 
-lll.Olg.e 


Sp 


fixed 


fixed 


0-015_o.oo5 


U.U1^_Q QQ4 


c 

OSp 


nxed 


fixed 


100 A +16.0 

-133.4I23.9 


1 OA 0+16.6 
-139.8_23.5 


Sv 


fixed 


fixed 


A A A 1 A AAC 

0.049 ± 0.005 


0.048 ± 0.005 


r 


nxed 


fixed 


/in /I +22.2 
-49.4_i86 


rr+21.5 
-47.7Ii8j 


A 

A 


A OA*? 1 A A1 O 

0.807 ± 0.013 


U.oU / zt U.Uio 


n finn -l n ni o 
u.ouy zt u.uiz 


A OAA 1 A A1 

0.809 ± 0.012 


P 


n 1 c;i _i_ n noc 


0.146 ±0.035 


0.116 ±0.030 


n inQ+0030 
u.iuy_g 


V 


0.401 ± 0.030 


0.400 ±0.030 


0.373 ± 0.029 


0.371 ±0.030 


xVdof 


20.7/8 


19.9/8 


44.6/30 


44.5/30 



Table 3: Fit results (l-o" ranges) of the theory parameters for Classes (A) and (B) in 
the two schemes defined in the text. The minimal value and the number of degrees of 
freedom (dof) are also given. The amplitudes are given in units of 10"^ eV, and the phases 
are in degrees. 

our fitting set from Class (A) to Class (B). They are correlated because they appear 
in combination in the physical amplitudes. 

The best fitted ratios between color-suppressed tree and tree amplitudes are 

(lA) (2A) (IB) (2B) 

Cv/Tv= 0.64 ±0.20 0.58 ±0.18 0.81 ±0.15 0.76 ±0.14 (3.1) 
Cp/Tp= 0.13 ±0.32 0.25 ±0.31 0.22 ±0.16 0.15 ±0.16. 

In the four schemes, the central values of the ratio Cp/Tp range from 0.13 to 0.25, 
agreeing with our naive expectation, even though one still cannot take them seriously 
due to the large errors coming from the uncertainty in Cp. On the other hand, the 
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central values for Cy jTy are significantly larger with less uncertainties. The value of 
Cy increases by about 40% from Set (A) to Set (B) though. The four schemes favor 
Cy jTy in the range of 0.58 ~ 0.76. As a comparison, the default parameter set of 
the QCDF approach |jl6[| gives 



CyjTy = 0.158 ± 0.109 and Cp/Tp = 0.20 ± 0.13, (3.2) 
The large Cy/Ty ratio is close to what we have found for C/T ~ 0.65 in the B 



decays to two pseudoscalars ^ ^ ^ ^ ^ H, Even though such values 



of C/T in the PP decays and Cy/Ty pose a challenge to perturbative calculations, 
they seem to follow the simple pattern of factorization in tree and color-suppressed 
tree amplitudes. 

The best fitted ratios between the QCD penguin amplitudes and the tree ampli- 
tudes are pretty stable among different schemes considered in this work. The ratios 
for the four schemes are given by 

(lA) (2A) (IB) (2B) 

Py/Ty= 0.06 ±0.01 0.06 ±0.01 0.06 ± 0.01 0.06 ± 0.01 (3.3) 
Pp/Tp= 0.12 ±0.01 0.12 ±0.01 0.11 ±0.01 0.11 ± 0.01 . 

In comparison, the ratio P/T ~ 0.21 in the PP modes The strong phase of Pp 
is the same as Tp within a few degrees, whereas that of Py is about 180° different. 



This agrees the expectation of Refs. ^ and reassures our previous finding 



2q] using old data. However, it is worth noting that in this work this solution is 



found even without invoking the B — > K*r) decays. The QCDF default values are 



le 



Py/Ty = 0.035 ± 0.017 and Pp/Tp = 0.032 ± 0.006 , (3.4) 

and also favors an opposite phase between Pp and Py amplitudes. Such a phase 
difference is due to the chiral enhancement that results in a sign flip in the effective 
coefficients for the QCD penguin amplitudes. Note, however, that the magnitudes of 
the QCD penguin amplitudes derived in QCDF are significantly smaller than what 
we find. It has been noticed that they cannot account for some large branching ratios 
in the QCD penguin-dominated modes ITH . 



The strong phase between Pp and Py is about 180°, with the former roughly in 
phase with Tp. Such a phase difference produces maximal constructive or destructive 
interference effects in decay modes that involve both of them. Since the relative 
phases among the tree- and penguin-type amplitudes are trivial {i.e., ~ 0° or 180°), 
as will be seen later, we generally do not expect large direct CP asymmetries in the 
decay modes involving only them. 

For the EW penguin amphtudes, the constraint on Pew,p is better than PEwy 
in Set (A). Nevertheless, the constraint on PEwy improves in Set (B). We note that 
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the strong phases of Pew,p and Pew,v are significantly different from those of Pp and 
Py, unhke the assumption made in Ref. It is interesting to notice that Pew,v 
increases by about 50% from fits of Set (A) to fits of Set (B). At the same time, the 
uncertainty in the strong phase associated with PEwy improves. In Set (B), Pswy 
is about 3 times larger than Pew,p- To one's surprise, PEwy is unexpectedly large, 
in line with Cy- 

As to the singlet penguin amplitudes, we find that Sp is about 3 times smaller 
than Sv- This partly justifies the ignorance of the former made in Ref. p8|, in 
view of the Okubo-Zweig-lizuka (OZl) rule. Moreover, if one compares the central 
values, the Sp amplitude has a strong phase in roughly the opposite direction of Pp 
and subtends a nontrivial angle from Cp. The Sy amplitude has a ~ 220° phase 
shift from Py and deviates from Cy by about 30°. It is interesting to note that the 
physical amplitude sp has a completely constructive interference between Sp and 
Pew,p/^- Also, both types of singlet penguin amphtudes are about half the sizes of 
the corresponding EW penguin amplitudes. 

Here we describe qualitatively how some of the theory parameters are fixed by 
data, thereby explaining their associated uncertainties. For this, we temporarily 
concentrate on the modes without involving singlet penguin amplitudes. But the 
argument can be easily extended to all modes. In our fits, the determination of Pp 
and Py is most precise because they can be directly extracted from the strangeness- 
changing K*^^^ and p^K^ modes. The next precise parameters are the 
magnitudes of tree amplitudes and their phase shifts relative to the QCD penguin 
amplitudes. They are fixed mainly by the strangeness-conserving p^-K^ and 
to some extent by the strangeness-changing B^ p~K^ and K*'^p^ modes. Since 
no direct CP asymmetry is observed in these modes, the relative strong phases are 
seen to be trivial. 

As the color-suppressed and EW penguin amplitudes of the same type (subscript 
P or V) always show up in pairs in the physical processes, the determination of 
their sizes and strong phases becomes trickier. This is because the color-suppressed 
amplitudes dominate in the AS* = processes, whereas the EW penguin amplitudes 
play more role in the \AS\ = 1 decays. This explains why Cy is better determined 
whereas PEwy is not, for B{B^ p+7r°) is more precise than B{K*'^n^). Likewise, 
the precision on Cp is worse than Pew,p because the combination of B{B^ p^K^) 
and B{B^ p^K^) is better than i3(p°7r+). 

Since the singlet penguin amplitudes are loop-mediated, they are better con- 
strained by the \^S\ = 1 decay modes. Currently, both charged and neutral (pK 
modes have consistent branching ratios and direct CP asymmetries. This basically 
fixes the magnitude and phase of Sp. In contrast, Sy is constrained in a more 
involved manner through interference with other amplitudes. 

We note in passing that in Class (A), we have also found other sets of parameters 
that render smaller Xmm ^^e fits. They are not listed in the tables because they 
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are not favored once the modes involving the singlet penguin amplitudes are taken 
into account. A distinctive feature of such solutions from the above-mentioned ones 
is that either the relative strong phase between Pp and Py is close to zero or that 
between Tp and Ty is close to 180°. In the former case, an interesting feature is 
that the ratios Cp/Tp = 0.57 ± 0.43 and Cv/Ty = 0.49 ± 0.12 in Scheme 2. They 
become comparable to each other, but still much larger than the usual perturbative 
expectation. In the latter case, we obtain a somewhat small p = 0.08. 

In Fig. m, we show the contours of the (p, f/) vertex at the 1-a and 95% confidence 
level (CL) obtained using Scheme 2. The left plot uses a fit to modes without 
involving the singlet penguin amplitudes. In this case, our favored region of the 
vertex is slightly higher than that given by the CKMfitter and UTfit [0. The 



right plot is a global fit to all the VP modes. Comparing to the left plot, we see 
that the favored region shifts lower and to the left on the p-f] plane. In this case, 
the preferred value of f3 agrees with other methods, while the value of 7 is slightly 
larger. The best fitted three angles in the UT are 



for Scheme (2A), and 



a = 


(83 ±8)° 


or 


72° 


<a< 99°(95%CL) 


/3 = 


(26 ±2)° 


or 


18° 


< /3 < 30°(95%CL) 


7 = 


(71 ±5)° 


or 


62° 


< 7 < 78°(95%CL) 


and 










a = 


(84 ±6)° 


or 


77° 


<a< 95°(95%CL) 


f3 = 


(23 ±2)° 


or 


18° 


<P< 23°(95%CL) 


7 = 


(73 ± 4)° 


or 


67° 


< 7 < 81°(95%CL) 



(3.5) 



(3.6) 
for Scheme (2B). 

The best-fitted UT vertex from the VP modes is highly consistent with the one 
from the PP modes. When the fits do not involve singlet penguin amplitudes, both 
VP and PP data favor a slightly larger 7 ^ 70° and a larger (3 ^ 26°. After including 
the modes involving the singlet penguin amplitudes, the best fitted 7 becomes even 
larger while /3 reduces to the value consistent with the B {cc)Ks measurements. 



4. Discussions 



There are two sets of decay modes that can provide a good test for the SU(3) sym- 
metry. One set contains the B+ K*'^n+, B+ K*^ K+ , B^ K*^ , and 
Bs K*^lf modes. The other set contains the B+ p+K^, B+ K*+lf , 
B'^ — s> K*^K , and Bs ^ K modes. They all involve only the Pp or Py ampli- 
tude, where we have neglected the Pewp ^ewv amplitude in the analysis as said 
before. However, this argument still applies if the color- suppressed EW penguin am- 
plitude is included because it scales in the same way as the QCD penguin amplitude. 



- 11 - 




Figure 1: The 1-a and 95% CL contours of the {p, ff) vertex obtained from a fit using the 
observed VP modes that do not involve the singlet penguin amplitudes (left) and using all 
of the observed VP modes (right), both assuming the exact flavor SU(3) symmetry. The 
l-fj range given by the CKMfitter is indicated by the cross. 



Currently, only the -B+ K*°7r+ and B 



P 



'K^ modes are observed, and their 



branching ratios are measured at (9(10"^) level. It is thus very helpful to measure 
any of the K*K modes in this respect. Using the fit results in Scheme (2 A) and in 
units of 10~^, we predict the branching ratios for the first set to be 10.64 ± 0.82, 
0.50 ± 0.05, 0.47 ± 0.05, and 9.11 ± 0.70, respectively. The branching ratios for the 
second set are 6.08 ± 0.79, 0.29 ± 0.04, 0.27 ± 0.04, and 5.21 ± 0.68 in units of 10"^ 
respectively. These B^^d K*K modes are somewhat difficult to measure due to the 
Cabibbo suppression. However, the B^ — > K*K modes should be within the reach 
of the LHCb and Tevatron Run-II experiments. 



Although the B^ 



!)7r+ and B^ 



i modes directly constrain the size 



of sp, their branching ratios are expected to be about (9(10^^) or smaller. Therefore, 
they are beyond the current probes. 

In the following, we would like to point out some persistent problems encountered 
in our fits to the current data. In Ref. [ESI, the rate difference relations [HHl: 



r(5° 
r(s° 



P TT 



p+vr-) 



r(5° 
r(5° 



P TT 



p-7r+) 



A [r(5° 

JK 



fp 



K* 



[r{B' 



^ K*-7i+) 



r(5'^ 



p-K 



(4.1) 
(4.2) 



have been found to be barely and loosely obeyed, respectively, by the data at that 
time. Using the current data and in terms of the branching ratios, Eqs. (O!) and 
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Mode 


BR (xlO ^) 


Acp 


s 




Bu,d p"7r+ 


16.59 ±4.01 (-0.09) 


-0.042 ± 0.041 (2.698) 


0.010 ± 0.173 (- 


-0.384) 


1 


7.52 ± 1.97 (0.05) 


0.049 ± 0.086 (-1.576) 


0.082 ± 0.166 (- 


-0.171) 


n n 


1.97 ± 0.94 (0.06) 


0.035 ± 0.179 (-) 


-0.064 ± 0.297 


[-) 




10.94 ±3.87 (-0.03) 


-0.011 ±0.193 (0.277) 


— 




p°7r+ 


8.81 ± 2.61 (-0.11) 


-0.121 ±0.090 (0.407) 


— 






0.47 ± 0.05 (-) 


(-) 


— 






0.27 ±0.04 (-) 


(-) 


— 






0.50 ± 0.05 (0.94) 


(-) 






K*+K^ 


0.29 ± 0.04 (-) 


(-) 








8.89 ± 1.13 (-0.29) 


0.094 ± 0.094 (0.926) 






p^fC" 


5.65 ± 1.21 (-0.26) 


0.076 ± 0.031 (-0.331) 


0.824 ± 0.047 (- 


-0.822) 


4- ry-n 


6.08 ± 0.79 (1.33) 


(-0.706) 






p^K"" 


3.80 ± 0.96 (0.03) 


0.382 ±0.119 (0.401) 


— 






6.59 ± 3.85 (-5.99) 


-0.330 ± 0.120 (1.500) 








8.87 ± 0.76 (1.30) 


-0.043 ±0.075 (-1.882) 








10.64 ± 0.82 (-0.80) 


(-0.339) 








7.00 ± 4.49 (-0.04) 


-0.081 ± 0.272 (0.418) 








6.89 ± 1.81 (-) 


0.049 ± 0.086 (-) 


— 




p°/C° 


0.39 ± 0.07 (-) 


0.929 ± 0.195 (-) 


-0.357 ± 0.528 


(-) 




15.22 ±3.68 (-) 


-0.042 ± 0.041 (-) 


— 






2.60 ± 1.25 (-) 


-0.134 ±0.328 (-) 


— 






7.45 ± 0.93 (-) 


0.085 ± 0.084 (-) 


— 






8.16 ± 0.70 (-) 


-0.041 ± 0.072 (-) 








5.21 ±0.68 (-) 


o(-) 


— 






9.11 ±0.70 (-) 


o(-) 






Table 4: Predicted B^^^s decay observables 


in Scheme (2A). Numbers in 


the parentheses 





are the pulls of theory predictions from the current experimental data. 

( [4.2| ) give in units of 10~^, respectively, 

-3.9 ± 2.0 = 2.1 ±0.9 , (4.3) 
2.1 ± 1.9 = -4.9 ±2.2 . (4.4) 

The first one is still not obeyed at about 2.7(t level. This difference comes from the 
CP asymmetries of p^ir^ and p^K^, both at about 2a level. To further check 

the equality in the second equation relies on more precise determinations in the CP 
asymmetries of B^ — > K*~^'n'~ and 5" —>■ K*^7t~ . 

Another problem is B{B^ — » p^t]') /B{B^ P~^v) — 1-3 ± 0.5, which is very 
different from our expectation of about 1/2 based upon the mixing angle we assume 
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for T] and rj' and assuming that sy is negligible for AS* = decays. The problem 
comes from the large branching ratio of P^v' i indicated by the pull in 

Table ^. A similar relation can be found for B{B^ p^ri)/B{B^ P^v')i ^{B^ 
u;°r))/B{B^ u^r]'), B{Bs p^r])/B{Bs p^r]'), and B{Bs Ljr])/B{Bs lot]') 
too. However, these modes may be difficult to measure. 

A new problem would occur between the B^ p^K^ and ujK^ modes that differ 
by \pls'p if s'p is vanishingly small. In that case, the ratio of their branching ratios 
should be close to 1 However, the current data comes down to 0.57±0.08. With 
the fitted Sp ~ 140 eV, the predicted ratio is ~ 0.61. Consequently, a non-vanishing 
Sp is preferred. 

Another puzzle comes from the CP asymmetry of B^ K*^-)] because it is 
measured at an almost Aa level. This is quite different from a closely related mode, 
B^ K*'^ri, whose CP asymmetry is consistent with zero. Their values should not 
be so different because they only differ by a small tree amplitude. 

We make predictions for the observables of all the B^, B^ and Bg decays using 
the extracted parameters given in Table ^. In Table ^ we only include modes without 
involving the singlet penguin amplitudes as they are based on Scheme (2A). Table ^ 
and Table ^ cover all the decay modes as they are based on Scheme (2B). The column 
of ^cp refers to either the direct CP asymmetry or A in Eq. (|2.5| ) of the corresponding 
mode. The numbers in the parentheses are calculated pulls of the theory predictions 
from experimental observations. They indicate the Ax^ contributions of individual 
quantities. 

Several observables in Table ^ have pulls larger than, say 1.5. Most of them are 
in the CP asymmetries. It is less clear about their importance as current precision 
on these data points is not satisfactory. We are then left with two branching ratio 
predictions with large pulls. The problem with p'^rj' has been mentioned above. As 
commented before, we do not include the branching ratio and CP asymmetry of the 
tt" in the fits of this work. Its predicted branching ratios in Tables ^ and 
^ based on the best fits are quite different from the current quotes of averages in 
Table and need further experimental confirmation. 

In Table |, our predictions of B{B^ p^r]) = 1.87 ± 0.64 and B{B° un^) = 
2.82 ± 0.99 are larger than the current upper bounds of 1.5 and 0.5, respectively, in 
units of 10~^. The branching ratio predictions of the other yet-measured modes are 
all consistent with current 95% upper bounds. 

For the B^ decays, we predict large direct CP asymmetries Acp{K rj) ^ 0.73 
and Acp{K rj') ^ —0.79, a result of interference between the large color- suppressed 
amplitude Cy and the QCD penguin amplitudes. We also predict large branching 
ratios, in unit of 10-^ B{<pi) - 8.47, B{K*~ti+) ~ 15.21, B{K*^K^) ~ 8, and 
B{K*^K ) ~ 9.54. In these modes, the branching ratios can reach (9(10^^) or more, 
as they involve either Ty for AS" = or Pp for \AS\ = 1 transitions. 
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5. Summary 



We have updated the global analysis of charmless B VP decays in the framework 
of flavor SU(3) symmetry using the latest experimental data. Moreover, we consider 
different SU(3) breaking schemes for the sizes of flavor amplitudes based upon factor- 
ization assumption. Our result shows that the symmetry-breaking scheme (Scheme 
2 defined in the text) is favored by the fits, but its difference from the exact sym- 
metry scheme (Scheme 1) is small. The UT vertex (p, f/) extracted using these modes 
is consistent with our previous analysis using the PP modes , and also agrees with 
other methods within errors []36| , However, we note that a slightly larger weak 
phase 7 is favored by our global analysis. 

In the fits to modes without involving the singlet penguin amplitudes, we note 
that there are two sets of solutions with minimal values. In one set, the Pp and 
Pv amplitudes have almost the same strong phases. In the other set, they have 
almost opposite strong phases. The latter is favored when one also includes modes 
involving the singlet penguin amplitudes. Moreover, we find in the latter case that 
the ratio Cy /Ty is about 0.6 - 0.7, similar to the C /T ratio in the PP modes. 
Correspondingly, the Pew,v and Sy amplitudes are unexpectedly large. These facts 
are seen to be a challenge to perturbative approaches. 

We point out that a set of decay modes that involve only the QCD penguin 
amplitude can be used to test our flavor SU(3) assumption. Among those modes, 
the Bg K*^K and K modes should be within the reach of the LHCb and 
Tevatron Run-II experiments. 

We also mention the persistent problems that the CP rate differences in B^ 
p^7r+ and in B^ p^K^ do not follow our expectation from factorization and that 
the observed branching ratio of B^ p^rj' is too large to be accommodated in 
our approach. Further investigations of B{B^ K*^ti^) and Acp{B^ — > K*^7]) are 
required. 

Based on our best fits, we calculate all observables in the B VP decays. The 
part for Bg decays is particularly useful because currently no such observables have 
been observed yet and our results serve as predictions to be compared with. 
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Mode 


BR (xlO-6) 




Acp 


S 






16.57 ±4.18 (-0.08) 


-0.038 ± 0.041 (2.630) 


0.070 ±0.166 (- 


-0.843) 




7.32 ± 1.98 


;o.2i) 


0.024 ± 0.072 (-1.363) 


0.084 ±0.160 (- 


-0.187) 




1.91 ±0.79 ( 


;o.i9) 


0.259 ± 0.148 (-) 


0.115 ± 0.249 (- 


-) 




11.12 ±2.99 


(-0.15) 


-0.026 ±0.128 (0.415) 


— 




n 1 


8.27 ±2.42 


;o.4i) 


-0.192 ± 0.099 (0.977) 


— 




n 


1.87 ±0.64 


:-) 


0.109 ± 0.153 (-) 


-0.336 ± 0.199 


(-) 




0.52 ±0.15 


:-) 


-0.396 ± 0.291 (-) 


-0.587 ± 0.222 


(-) 




7.16 ±2.03 


-0.26) 


0.165 ±0.103 (-0.502) 


— 






3.79 ±0.98 


1.63) 


-0.071 ± 0.240 (0.110) 


— 






2.82 ± 0.99 


;-) 


0.293 ± 0.132 (-) 


-0.094 ± 0.216 


(-) 




7.02 ± 2.23 ( 


-0.25) 


0.020 ± 0.075 (-0.993) 


— 




CUT) 


1.27 ±0.51 ( 


;-) 


-0.016 ±0.179 (-) 


-0.360 ± 0.227 


(-) 


ijjrl 


0.76 ±0.25 


;-) 


-0.624 ±0.285 (-) 


-0.511 ± 0.302 


(-) 




0.02 ±0.01 


-) 


(-) 


(-) 






0.04 ±0.02 


-) 


(-) 


— 






0.01 ±0.01 


-) 


O(-) 


(-) 




H 


0.01 ±0.00 


-) 


O(-) 


O(-) 






0.52 ± 0.05 


;-) 


O(-) 


— 




K*°K° 


0.31 ± 0.04 


;-) 


O(-) 


— 




K*^K+ 


0.55 ± 0.05 ( 


;o.67) 


(-) 


— 






0.33 ± 0.04 ( 


;-) 


O(-) 


— 






9.21 ± 1.04 


;-o.6i) 


0.082 ± 0.089 (1.128) 






p^K^ 


5.06 ±1.10 


;o.36) 


-0.041 ± 0.045 (0.072) 


0.766 ±0.052 (- 


-0.598) 


p-^K^ 


6.70 ±0.74 


;o.9o) 


(-0.706) 


— 




p^K+ 


4.02 ±0.82 


;-0.44) 


0.382 ±0.126 (0.398) 


— 






4.62 ± 1.01 


;o.63) 


0.033 ± 0.048 (1.690) 


0.700 ±0.054 (- 


-1.040) 


T J- 1 


6.64 ±1.27 


;o.i3) 


0.029 ± 0.092 (-0.190) 








7.43 ± 1.21 


;o.79) 


(1.533) 


0.737 ±0.043 (- 


-1.699) 




7.96 ±1.30 


;o.53) 


(0.773) 


— 






13.85 ±4.76 


(-16.36) 


-0.294 ± 0.078 (1.201) 


— 






9.57 ±0.72 


;o.66) 


-0.019 ±0.057 (-2.104) 


— 






11.14 ±0.77 


(-1.43) 


(-0.339) 








7.09 ± 3.11 


;-0.08) 


-0.151 ± 0.164 (0.660) 








16.72 ± 2.44 


(-0.82) 


0.162 ± 0.049 (0.560) 








4.16 ± 1.56 ( 


-0.30) 


0.159 ±0.150 (-0.954) 








17.30 ±2.58 


(1-25) 


0.070 ± 0.064 (-0.837) 






K*+r]' 


4.34 ±1.64 


;o.28) 


-0.027 ±0.228 (0.933) 







Table 5: Predicted d decay observables in Scheme (2B). Numbers in the parentheses 
are the pulls of theory predictions from the current experimental data. 
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Mode 


BR (xlO-6) 




Acp 




s 


p^rj 


0.21 ±0.14 ( 


-) 


-0.156 ±0.123 (- 


) 


-0.731 ±0.092 (-) 


p^'i 


0.42 ± 0.26 ( 


-) 


-0.156 ±0.123 (- 


) 


-0.731 ±0.092 (-) 


p-K+ 


6.71 ± 1.81 ( 


-) 


0.024 ±0.072 (-) 




— 


p'^K^ 


0.24 ±0.10 ( 


-) 


-0.128 ±0.773 (- 


) 


0.926 ±0.283 (-) 


LOT] 


0.07 ±0.06 ( 


-) 


0.243 ±0.234 (-) 




-0.624 ±0.195 (-) 


UJT]' 


0.13 ±0.12 ( 


-) 


0.243 ±0.234 (-) 




-0.624 ±0.195 (-) 




0.27 ±0.14 ( 


-) 


0.302 ±0.629 (-) 




-0.856 ±0.331 (-) 




2.80 ± 1.80 ( 




-0.250 ±0.121 (- 


) 


-0.451 ±0.131 (-) 


077 


2.35 ± 1.53 ( 




-0.073 ±0.142 (- 


) 


-0.341 ±0.174 (-) 


077' 


8.47 ±2.55 ( 




0.096 ±0.061 (-) 




-0.626 ±0.054 (-) 




0.44 ± 0.07 ( 




OH 




o(-) 


K*-7r+ 


15.21 ±3.83 


(-) 


-0.038 ±0.041 (- 


■) 




kV 


4.27 ± 1.36 ( 




-0.064 ±0.146 (- 


■) 






3.26 ±0.93 ( 




0.730 ±0.108 (-) 








1.99 ±0.47 ( 




-0.794 ±0.191 (- 


■) 




K*-K+ 


7.79 ± 0.86 ( 




0.073 ±0.079 (-) 






K*+K- 


8.79 ±0.66 ( 




-0.018 ±0.054 (- 


) 






5.74 ±0.63 ( 




OH 








9.54 ±0.66 ( 




OH 







Table 6: Predicted Bg deeay obscrvables in Seheme (2B). Numbers in the parentheses are 
the puUs of theory predictions from the current experimental data. 
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